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Abstract: Consider (for simplicity) two one-dimensional semi-infinite leads 
coupled to a quantum well via time dependent point interactions. In the re- 

Smote past the system is decoupled, and each of its components is at thermal 
equilibrium. In the remote future the system is fully coupled. We define and 
compute the non equilibrium steady state (NESS) generated by this evolution. 
£S) . We show that when restricted to the subspace of absolute continuity of the fully 

J> ' coupled system, the state does not depend at all on the switching. Moreover, we 

show that the stationary charge current has the same invariant property, and 
C*") ■ derive the Landau-Lifschitz and Landauer-Buttiker formulas. 

ON ■ 

cn ; 

00 ' 1- Introduction 

o: 

t^- , The goal of this paper is to construct and study non equilibrium steady states 

for systems containing quantum wells, and to describe the quantum transport 
of electrons through them. Even though our results can be generalized to higher 
•i"H , dimensions, we choose for the moment to work in a (quasi) one dimensional 

' setting; let us describe it in some more detail. 

A quantum well consists of potential barriers which are supposed to confine 
particles. On both sides of the barriers are reservoirs of electrons. Carriers can 
pass through the barriers by tunneling. We are interested in the carrier trans- 
port through the barriers, as well as in the carrier distribution between these 
barriers. Models of such type are very often used to describe processes going on 
in nanoelectronic devices: quantum well lasers, resonant tunneling diodes, and 
nanotransistors, see [28] , 

The quasi one dimensional geometry assumes that the carriers can freely move 
in the plane orthogonal to the transport axis, but these degrees of freedom are 
integrated out. Thus we are dealing with an essentially one-dimensional physical 
system. To describe such a system we consider the transport model of a single 



2 



H.D. Cornean, H. Neidhardt, V. A. Zagrcbnov 



band in a given spatially varying potential v, under the assumption that v and 
all other possible parameters of the model are constant outside a fixed interval 
(a, b), see [HUMID . 

More precisely, in the Hilbert space 9) := L 2 (M.) we consider the Schrodinger 
operator 

^^■■=-\iik)i s{x)+v{xmxl xeR ' M 



with domain 



Dom (H) := {/ e W 1 ' 2 ^) : -^f 6 W l,2 (R)}. (1.2) 



It is assumed that the effective mass M(x) and the real potential V(x) admit 
decompositions of the form 

{m a x £ (— oo, a] 
m(x) x€(a,b) , (1.3) 
mij x £ [6, oo) 

< m a , To;, < oo, m(x) > 0, x £ (a, 6), m + m 1 (b) £ L°°((a, b)), and 



w a a; £ (— oo, a] 

a; £ (a,b) , w a > u b , (1.4) 
Vb x € [b, oo) 



v a , Vb £ M, w £ L°°((a, fe)). The quantum well is identified with the interval (a, 6), 
(or physically, with the three-dimensional region (a, b) xR 2 ). The regions (— oo, a) 
and (6, oo) (or physically (— oo, a) x R 2 and (b, oo) x M 2 ), are the reservoirs. 

Schrodinger operators with step-like potentials were firstly considered by Bus- 
laev and Fomin in [8]. For that reason we call them Buslaev-Fomin operators. 

The inverse scattering problem for such Buslaev-Fomin operators was subse- 
quently investigated in [Ilf2l[3TfT0lfTTl [T8] . 

In order to rigorously describe quantum transport in mesoscopic systems, 
these operators were firstly used by Potz, see [29j . In [6], the Buslaev-Fomin 
operator was an important ingredient for a self-consistent quantum transmitting 
Schrodinger-Poisson system, which was used to describe quantum transport in 
tunneling diodes. In a further step, this was extended to a so-called hybrid model 
which consists of a classical drift-diffusion part and a quantum transmitting 
Schrodinger-Poisson part, see [7]. Hybrid models are effective tools of describing 
and calculating nanostructures like tunneling diodes, see [5]. 

To obtain a self-consistent description of carrier transport through quantum 
wells, one needs to know the carrier distribution between the barriers in order to 
put it into the Poisson equation for determining the electric field. Semiconductor 
devices are often modeled in this manner, see [1 7i[23][3D, . Important for that is 
a relation which assigns to each real potential v £ L°°((a, b)) a carrier density 
u £ L 1 ((a, 6)). The (nonlinear) operator doing this is called the carrier density 
operator and is denoted by 

TV(-) : L°°((o, b)) — » ^((a, 6)), 7V» = u. 
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The problem of defining carrier density operators is reduced to the problem 
of finding appropriate density operators g. 

Definition 11 A bounded non-negative operator g in L 2 (R) is called a density 
operator or a state if the product gM(x(a,b)) * s a trace class operator, where 
M(x(a,b)) * s ^ e rnultiplication operator induced in L 2 (R) by the characteristic 
function X{a.b) of the interval (a, b). 

We note that in general a non-negative bounded operator is called a state if 
the operator itself is a trace class operator and is normalized to one, that is, 
Tr(fj) = 1. In our case these conditions are relaxed to the condition that the 
product gM{x( a .b)) has to be trace class. 

This weakening is necessary since we are interested in so-called steady density 
operators or steady states for Hamiltonians with continuous spectrum. 

Definition 12 A state g is called a steady state for H if g commutes with 
H, i.e. g belongs to the commutant of the algebra generated by the functional 
calculus associated to H . A steady state is an equilibrium state if it belongs to 
the bicommutant of this algebra. 

Thus if H admits continuous spectrum, then a steady state cannot be of trace 
class unless it equals zero on the subspace of absolute continuity. 

To give a description of all possible steady states, one has to introduce the 
spectral representation of H. Taking into account results of [Bj, it turns out that 
the operator H is unitarily equivalent to the multiplication M induced by the 
independent variable A in the direct integral L 2 (R, f)(A), v), 



C, A e (-co, v a ] 
C 2 , Xe(v a ,oo) 



f)W:=Ca .A. 17 J , (1-5) 



and (with the usual abuse of notation) 

N 
3=1 

where it is assumed v a > u&, and {Ajj^j denote the finite number of simple 
eigenvalues of H which are all situated below the threshold Vb- We note that 

L 2 (R, f)(A), v) ~ ©f =1 C © L 2 ([v b , t)„],C) © L 2 ((v a , oo), C 2 ). 

The unitary operator <P : L 2 (R) — ► L 2 (R, f)(A), v) establishing the unitary 
equivalence of H and M is called the generalized Fourier transform. 
If g is a steady state for H, then there exists a z^-measurable function 

R 3 A h-> p(X) S B(t){\)) 

of non-negative bounded operators in t)(A) such that v — sup AgR |[p(A) |[«B(t)(A)) < 
oo and g is unitarily equivalent to the multiplication operator M(p) induced by 
p via the generalized Fourier transform 

g = <P~ 1 M{p)$. (1.6) 
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The measurable family {/5(A)}asr is uniquely determined by the steady state g 
up to a v-zero set and is called the distribution function of the steady state. In 
other words, there is an one-to-one correspondence between the set of steady 
states and the set of distribution functions. When p is an equilibrium state, then 
/5(A) must be proportional with the identity operator in f)(A), hence p must be 
a function of H. Let us note that the same distribution function can produce 
quite different steady states in L 2 (R) . This is due to the fact that the generalized 
Fourier transform strongly dependents on H, in particular, on the potential v. 

Having a steady state g for H one defines the carrier density in accordance 
with |6J as the Radon-Nikodym derivative of the Lebesgue continuous measure 
E(w) 

E(w) := Tr( Q M( Xu )) 

where a; is a Borel subset of (a,b). The quantity E(w) can be regarded as the 
expectation value that the carriers are contained in lo. Therefore the carrier 
density u is defined by 

dx dx 

The carrier density operator Af p (-) : L°°((a, b)) — > L}{{a, b)) is now defined as 

N p {v) := u p {x) (1.7) 

where v £ L°°((a,b)) is the potential of the operator H. The steady state g is 
given by 1(1.60 . 

Therefore the self-consistent description of the carrier transport through 
quantum wells is obtained if there is a way to determine physically relevant 
distribution functions p. One goal of this paper is to propose a time-dependent 
procedure allowing to determine those functions. 



1.1. The strategy. Let us describe the strategy. We start with a completely de- 
coupled system which consists of three subsystems living in the Hilbert spaces 

fi a :=L 2 ((-oo,a]), fix := L 2 (l), fi b := L 2 ([b,<x)) (1.8) 

where I = (a, b). We note that 

ft = fia ® fix ® fib- (1.9) 
With Sj a we associate the Hamiltonian H a 

( Ha f)(x):=--^-^f(x)+v a f(x), (1.10) 
2m a dx z 

I e Bom(H a ) := {/ £ W^ 2 ((-ao,a)) : /(a) = 0} (1.11) 
with fix the Hamiltonian Hi, 

^^ : =-^^)l^ + ^)' ^ 



f E Dom (2fc) := {/ e ^-(X) : ^ = /(6) = } 



(1.13) 
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and with 9j b the Hamiltonian H b 



(^/)(*):=-^/(*)+W(*). (1-14) 
/ e Dom(H b ) := {/ e W 2 ' 2 {{b^) : f(b) = 0}. (1.15) 



In f) we set 



H D :=H a ®H x ®H b (1.16) 

where the sub-index "D" indicates Dirichlct boundary conditions. The quantum 
subsystems {$j a , H a } and {$jb, H b } are called left- and right-hand reservoirs. The 
middle system {Sjx, Hx} is identified with a closed quantum well. We assume 
that all three subsystems are at thermal equilibrium; according to Definition 
I12[ the corresponding sub-states must be functions of their corresponding sub- 
Hamiltonians. The total state is the direct sum of the three sub-states. 

One example borrowed from the physical literature, which takes into account 
the quasi one dimensional features of our problem is as follows. Assume the same 
temperature T. The equilibrium sub-states are g a , qx and g b where: 

Qa'=fa(H a -fl a ), Qx := fx(Hj - fa), Qb ■= %{H b - fl h ) (1-17) 

where 

f«(A) := c a ln(l + e-? x ), f x (A) := c x ln(l + e"^), f b (A) := c b ln(l + e^ A ), 

A G R, (3 := l/(fcT), k is the Boltzmann constant, fj, a and fi b are the chemical 
potentials of left- and right-hand reservoirs and fix the chemical potential of the 
quantum well. The constants c a , cx and c b are given by 

r — qm * a r — qm * X r — 9 ^ MK\ 

7rp 7rp it p 

where m*, and mj* are the electronic effective masses appearing after inte- 
grating out the orthogonal degrees of freedom (see for more details [151115] ). 
We set 

QD ■= Qa © Qx® Qb- (1.19) 

For the whole system {Sj,Hu} the state qd is a steady state because com- 
mutes with Ho (see Definition ! 12|) . In general, the state qd cannot be represented 
as a function of Hd which is characteristic for equilibrium states, but it is the 
direct sum of equilibrium sub-states. In any case, qd is a special non-equilibrium 
steady state (NESS) for the system {Sj, Hd}- Now here comes the main question: 
can we construct a NESS for {Sj, H} starting from qd? 

Let us assume that at t = — oo the quantum system {fj,ii/£>} is described 
by the NESS qd- Then we connect in a time dependent manner the left- and 
right-hand reservoirs to the closed quantum well {f)i, Hx}- We assume that the 
connection process is described by the time-dependent Hamiltonian 

H a (t) :=H + e- at S(x~a) + e- at 5(x-b), t e R, a > 0. (1.20) 

The operator H a (t) is defined by 

(H a (t))f)(x):=-~^--^-f(x) + V(x)f(x), f e Bom (H a (t)), (1.21) 
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where the domain Dom (H a {t)) is given by 

Dom (#<*(*)) := (1-22) 
jjf G W^(R) 
f e W^iR) : (^f')(a + 0) - (^/'Xo - 0) = e-"\f(a) 
(4jf')(b + 0)-(4jn(b-0) = e-^f(b) 

After a rather standard analysis, one can prove the following convergence in the 
norm resolvent sense: 

n- lim (HJt) - z)- 1 = (H D - z)- 1 (1.23) 

t — > — OO 

and 

n - lim (HJt) - z)' 1 = (H - z) -1 , (1.24) 

z G C \ R. Then we consider the quantum Liouville equation (details about the 
various topologies will follow later): 

*J^a(*) = ea(*)], * e R, (1-25) 

for a fixed a > satisfying the initial condition 

s- lim g a (t) = q d . 

t — * — oo 

Having found a solution g a (t) we are interested in the ergodic limit 

T 



i r 

Qa. = Jun 7f / Qa{t)dt. 



If we can verify that the limit g a exists and commutes with H, then g a is 
regarded as the desired NESS of the fully coupled system {$),H}. Inserting g a 
into the definition of the carrier density operator Af Pa we complete the definition 
of the carrier density operator. Finally, the steady state g a allows to determine 
the corresponding distribution function {/5 Q (A)}AeK- 



1.2. Outline of results. The precise formulation of our main result can be found 
in Theorem [35j and here we only describe its main features in words. 
We need to introduce the incoming wave operator 

W- := s- lim e im e- UHD P ac (H D ) (1.26) 

t — > — oo 

where P ac (Ho) is the projection on the absolutely continuous subspace S) ac (Hn) 
of H D . We note that h ac (H D ) = L 2 ((-oo, a]) © L 2 {[b, oo)). The wave operator 
exists and is complete, that is, W_ is an isometric operator acting from S} ac (Ho) 
onto f) ac (H) where Sj ac (H) is the absolutely continuous subspace of H (the range 
of P ac {H)). 

One not so surprising result, is that g a exists for all a > 0. In fact, if we 
restrict ourselves to the subspace f) ac (H), then we do not need to take the 
ergodic limit, since the usual strong limit exist. The surprising fact is that 

s- lim g a (t)P ac (H) = g a P ac (H) = W- PD WlP ac {H), (1.27) 

t—*OC 
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which is independent of a. 

The only a dependence can be found in g a P d (H), where P d (H) is the pro- 
jection on the subspace generated by the discrete eigenfunctions of H . But this 
part does not contribute to the stationary current as can be seen in Section 3J 
Here the ergodic limit is essential, because it kills off the oscillations produced 
by the interference between different eigenfunctions. 

Note that the case a = oo would describe the situation in which the coupling 
is suddenly made at t = and then the system evolves freely with the dynamics 
generated by H (see [TH] , and the end of Section [3]) . 

The case a \ would correspond to the adiabatic limit. Inspired by the 
physical literature which seems to claim that the adiabatic limit would take care 
of the above mentioned oscillations, we conjecture the following result for the 
transient current: 

Conjecture 1. 

lim limsup\Tr{g a (t) P d (H)[H,x]}\ = 0, 

where x lS anv smoothed out characteristic function of one of the reservoirs. 

Before ending this introduction, let us comment on some other physical as- 
pects related to quantum transport problems. Many physics papers are dealing 
with transient currents and not only with the steady ones. More precisely, they 
investigate non-stationary electronic transport in noninteracting nanostructures 
driven by a finite bias and time-dependent signals applied at their contacts to 
the leads, while they allow the carriers to self-interact inside the quantum well 
(see for example [23], [25] and references therein). 

An interesting open problem is to study the existence of NESS in the Cini 
(partition- free) approach [9], [13], [14], [12]. Some nice results which are in the 
same spirit with ours have already been obtained in the physical literature |31j , 
[32], even for systems which allow local self- interactions. 

Now let us describe the organization of our paper. 

Section [2] introduces all the necessary notation and presents an explicit de- 
scription of a spectral representation of Ho and H . 

Section [3] deals with the quantum Liouville equation, and contains the proof 
of our main result in Theorem 1351 

In Section |4] we define the stationary current and derive the Landau-Lifschitz 
and Landauer-Biittiker formulas. 



2. Technical preliminaries 

2.1. The uncoupled system. Let us start by describing the uncoupled system, 
and begin with the left reservoir. The spectrum of H a is absolutely continu- 
ous and a(H a ) = o- ac (H a ) = [v a ,co). The operator is simple. The generalized 
eigenfunctions ip a (', A), A £ [v a , oo), of H a are given by 

, i u sin(2m a g a (A)(a; - a))) 

ip a (x, A) := , xe{-oo,a\, Xe[v a ,oo) 
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where 

Qa (A) 



2m n 



The system of eigenfunctions {ip a {'-> ^)}\e[v a .oc) is orthonormal, that is, one has 
in distributional sense 



dx %(x,X)ip a (x, fi) = 8(X- fx), X,/j,e[v a ,oo). (2.1) 

With the generalized eigenfunctions one associates the generalized Fourier trans- 
form ip a : L 2 ((-oo, a]) — > L 2 ([v a , oo)) given by 



(V a f)(X)= f(x)%(x,X)dx = f(x)xp a (x, X)dx 



f G L 2 ((— oo,a\). Using (|2.1|) a straightforward computation shows that the 
generalized Fourier is an isometry acting from L 2 ((— oo, a]) onto L 2 ([v a , oo)). The 
inverse operator if'" 1 : L 2 ([v a , oo) — > L 2 (— oo, a]) admits the representation 

(^/)(A)= f^(,,A)/(A)dA= rBin(2^ ga (A)(x- a) )) /(A)dA> 

f G L 2 ([v a ,oo)). Since -0a(-,A) are generalized eigenfunctions of H a one easily 
verifies that 

where M a is the multiplication operator induced by the independent variable A 
in L 2 ([v a , oo)) and defined by 

(M a f)(X) = A/(A), 
/ G Dom (Ma) := {/ G L 2 (K,oo)) : A/(A) G L 2 (K,oo))}. 

This shows that {L 2 ([v a , oo)), M a } is a spectral representation of H a . For the 
equilibrium sub-state g a = f a (H a — /i a ) one has the representation 

Qa = *a X M{U{- ~ Ha))*a 

where M(f a (- — /i a )) is the multiplication operator induced by the function f a (- — 

Ma)- 

Let us continue with the closed quantum well. The operator Hj has purely 
discrete point spectrum {£fc}fcgN with an accumulation point at +00. The eigen- 
values are simple. The density operator qj = \i(Hx — fJ-x) is trace class. One 
easily verifies that there is an isometric map ^x ■ L 2 ((a,b)) — > L 2 (K,C, vz), 
du%(X) = Y^k=i ^(A — £,k)dX, such that {L 2 (K, C, vx), Mx} becomes a spec- 
tral representation of Hx where Mj denotes the multiplication operator in 
L 2 (R,C,^x). 
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Finally, the right-hand reservoir. The spectrum of H b is absolutely continuous 
and a(H b ) = a ac (H b ) = [v b ,oo). The operator H b is simple. The generalized 
eigenfunctions ip b (-, A), A £ [v b , oo) are given by 

sm(2m b q b (\)(x - b)) 
ip b {x,X) = , xe[b,oo), Xe[v b ,oo) 

where 



The generalized eigenfunctions {V'bC'i A)}A6[i> b ,ao) perform an orthonormal sys- 
tem and define a generalized Fourier transform W b : L 2 ([b, oo)) — ► L 2 ([d{,,oo)) 

by 



(%f)(X):= f{x)^ b {x 1 X)dx = f(x)Mx,\)dx, 

J b Jb 

f G L 2 ([6, oo)). The inverse Fourier transform g^ -1 : L 2 ([w fc ,oo)) — > L 2 ([6, oo)) 
admits the representation 

/)W = / ih{x,X)f{X)dX= / /(A)dA, 



/ G L 2 ([v&, oo)). Denoting by M& the multiplication operator induced by the 
independent variable A in L 2 ([v b , oo)) we get 



/- 1 or H b =#- 1 



which shows that {L 2 ([v b , oo)), M b } is a spectral representation of H b . The equi- 
librium sub-state g b = f b (H b — fj, b ) is unitarily equivalent to the multiplication 
operator M(f b (- — fi b )) induced by the function F(- — (i b ) in L 2 ([v b , oo)), that is, 

g b =W b - 1 M(f b (--fx b ))W b . 



2.2. Spectral representation of the decoupled system. A straightforward compu- 
tation shows that the direct sum ^ = 9 a ®'i r i®^ / b defines an isometric map acting 
from L 2 (R) onto i 2 (R, f)(A), v D (X)), dv D {X) = £~ x f(A-&)dA + X[v»,oo)(A)dA> 
such that Ho becomes unitarily equivalent to the multiplication operator Mb 
defined in L 2 (R, f)(A), vd(X)) (see (|1.5|1 ). Here we slightly change the definition 
of f)(A) such that it re-becomes C when A hits an eigenvalue. This does not affect 
the absolutely continuous part. 

Hence {L 2 (R, f)(A), z/£>(A)), Md} is a spectral representation of Hp- Under 
the map \P the absolutely continuous part H D C — H a © H b of Hr> is unitar- 
ily equivalent to the multiplication operator M in L 2 (R, 1)(A), v D c }, dv D c (X) = 
X[v b ,oa){^)dX. Therefore {L 2 (M, f)(A), v D c ), M} is a spectral representation of H^f. 
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With respect to the spectral representation {L 2 (M., rj(A), I'd), M} the distri- 
bution function {p£>(A)} of the steady state go is given by 



'0, 

fx (A - 11%), 
p D (X) := < fb(A - fib), 



ffe(A-M6) 

, fa(A-/i«), 



A eR\a(H D ) 

A e a p (H D ) = a(H x ) 

A G [v b , v a ) \ a{H x ) 

A e [u a ,oo) 



We note that M(p D ) = ^g D ^~ 



2.3. The fully coupled system. The Hamiltonian H in (jl.ip was investigated in 
detail in [6]. If v a > «6, then it turns out that the operator H has a finite simple 
point spectrum on (— oo, Vb), on [vb,v a ) the spectrum is absolutely continuous 
and simple, and on [v Q ,oo) the spectrum is also absolutely continuous with 
multiplicity two. 

Denoting by {X p }p =1 the eigenvalues on (— oo, Vb), we have a corresponding 
finite sequence of L 2 -eigenfunctions {i(j(x, )}^ =1 . 

Moreover, one can construct a set of generalized eigenfunctions (f> a (x,X), 
x S R, A S [u a ,oo), and </>b(a;,A), x 6 M, A 6 [v&,oo) of _ff such that 
WiM)}Ae[« sA ) and {<fo(-, A), a (-, A)} Ae[ „ ai0o) generate a complete orthonor- 
mal systems of generalized eigenfunctions. More precisely: 



4> a {x, \)(f> a (x, n)dx = 5(\ - ji), A, ^ € [w Q , oo) 



A)^b(a;, ^)da; = 5(A — /i), A,^i £ [ut, oo) 



^ a (a;, \)</> b (x, n)dx — 0, A,^e[w a ,oo), 

see [6]. The existence of generalized eigenfunctions is shown by constructing 
solutions 4> a (x, A) and 4>b{x, A) of the ordinary differential equation 

i £ M, A € oo), p — a,b, obeying 



(j> a {x, A) = 
A G [ti a , oo), and 
<j)b(x,\) = 



e i2 ma q a (\)(x-a) + Saa ( X y- t 2rn a q a (\)(x-a) x £ fl ] 

Sba{X)e t2mbqbW{x ~ h] X € [6, oo), 



5 , a6(A)e- l2m " 9 "( A )( ;E - Q ) a; E (-co, a] 

e ~i2 mb q b {\)(x-b) + Sbb ^ e i2m, b q b {\)(x-b) % £ ^ 
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A G [vb,oo). The coefficients S aa (X) and Sbb(X) are called reflection coefficients 
while Sba{X) and S a b(X) are called transmission coefficients. The solutions </> a (A) 
and 4>b(X) define the normalized generalized cigcnfmictions of H by 

(f) b (x, A) := - — ^— — 4> b {x, A), x G M, AG[v b ,oc), 
47rq & (A) 

^o(ac,A):=- — -7 -ty 4> a {x, A), x G E, Ag[w q ,oo). 
47rg a (A) 

Having the existence of the generalized eigenfunctions one introduces the gener- 
alized Fourier transform <P : L 2 (R) — > i 2 (l)(K, f)(A), v) by 



(*/)(A) := / f(x)4>(x,\)dx, feL 2 (R), A € tr(tf), (2-2) 

where 



<f>(x,X) := 



<l>(x,\j) A G <7p(H), xel 
<p b {x, A) A G [vb,v a ), x G 
f <£i>(z, A) 



(2.3) 



_ \ o (x, A) 



A G [vb, oo), xeK, 



see [6]. The inverse generalized Fourier transform <P 1 : L 2 (R, f)(A), i/) — > L 2 
is given by 



($- 1 9)= / (g(A)^(x,A))„ (A) ^(A), (2.4) 

i 6 I, j 6 L 2 (M, f)(A), ^) where (•, is the scalar product in f)(A). Since ^ is 
an isometry action from L 2 (R) onto L 2 (M, f)(A), f) and 

M = 'PH'p- 1 

holds where M is multiplication operator induced by the independent variable 
in L 2 (R, t)(A), v) one gets that {L 2 (M, f)(A), ^), M} is a spectral representation 
of H. 

Under the absolutely continuous part H ac becomes unitarily equivalent 
to the multiplication operator M in L 2 (R, f)(A), v ac ), du ac — X[v b ,oo)dX. Hence 
{L 2 (R, f)(A), v ac ), M} is a spectral representation of i/ ac , as it was for H^f. 



2.4- The incoming wave operator. We have already mentioned that W- as de- 
fined in (|1.26p exists and is complete [33] . We will need in Section[4]the expression 
of the "rotated" wave operator $W—& which acts from L 2 (R, f)(A), v ac ) onto 
itself. By direct (but tedious) computations one can show that W- := <£W_<?' _1 
acts as a multiplication operator, which means that there is a family {W_ (A)},\£R 
of isometries acting from E)(A) onto f)(A) such that 



(W_/)(A) = W_(A)/(A), /GL 2 (R,f)(A),^ c ). 
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The family {W_(A)}asr is called the incoming wave matrix and can be explicitly 
calculated. One gets 

A g [v b ,v a ] 

(2.5) 




A G (V a , OO). 

Note that another possible approach to the spectral problem (and completely 
different) would be to construct generalized eigenfunctions for H out of those of 
Ho by using the unitarity of W— between their subspaces of absolute continuity 

and the formal intertwining identity (f> p (-, A) := W-ij} p (-, A). In this case W-(X) 
would always equal the identity matrix. 



3. The quantum Liouville equation 

The time dependent operators H a {t) from (|1 .20|) are defined by the sesquilinear 
forms f) a [t](; •), 

i) a [t](f,g)= (3.1) 
/ \f(x)J{x) + V(x)f(x)gjxj\ dx + e^'/WjW + e _at /(6)ffW, 

f,g G Dom ({)„[*]) := W 1,2 (R), t G K. Obviously, we have + r > I, 

T '■= \\V\\l°°(m.) + 1- F° r each t £ 1 the operator H a (t) can be regarded as a 
bounded operator acting from W ' (R) into W /_1,2 (R). Classical Sobolev embed- 
ding results ensure that (H a (t) + t)~ 1//2 maps L 2 (M.) into continuous functions, 
and it has an integral kernel G(x, x'; r) with the property that G(-, x'\ r) G L 2 (R) 
for every fixed x'. 

Let us introduce the operators B a : L 2 (R) — ► C and B b : L 2 (M) — > C 
defined by: 

(B a f):=[(H a (t)+r)- 1/2 f](a)= (3.2) 

/ G(a, x] r)f(x)dx, {B* a c)(x) = G(x, a; r), 
Js. 

and similarly for B b . The operators B* a B a and B b B b are bounded in L 2 (R) and 
correspond to the sesquilinear forms 

b a [t](f,g) := ((H a (t) + r)- 1 / 2 f)(a)((H a (t) + r)-^g)(a), 
f,g G Dom(baM) = £ 2 ( R ), and 

b b [t](f,g) := ((H a {t) + r)- 1 / 2 f)(b)((H a (t) + r)-^g)(b), 

f,g G Dom(bf,[t]) = L 2 (R), respectively. We define the rank two operator 

B:=B* a B a + B* b B b = (3.3) 
G(-, a; r)G(a, •; r) + G(-, 6; r)G(6, •; r) : L 2 (R) — > L 2 (R). 

The resolvent (H a (t) + t) _1 admits the representation 
(ff^il + r)- 1 = (H + T)- 1 ' 2 {I + e- at B)- 1 {H + r)- 1 ' 2 , t G R, a > 0. (3.4) 
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3.1. The unitary evolution. Let us consider a weakly differentiable map t3fn 
u(t) 6 W 1,2 (M.). We are interested in the evolution equation 

d 

i—u(t) = H a (t)u(t), iel, a>0. (3.5) 

where H a (t) is regarded as a bounded operator acting from W 1,2 (R) into 

By Theorem 6.1 of [26] with evolution equation (13 . 5|) one can associate a 
unique unitary solution operator or propagator {U(t, s)}(t iS )£r x r leaving in- 
variant the Hilbert space W 1 ' 2 (R). By Theorem 8.1 of [20J we find that for 
x,y £ W 1,2 (M.) the sesquilinear form (U(t, s)x,y) is continuously differentiable 
with respect t£l and s£l such that 

^- t (U(t,s)x,y) = -i(H a (t)U(t,s)x,y), x,y e W^fR), (3.6) 
^-(U(t,s)x,y)=i(H a (s)x,U(s,t)y), x,yeW^(I). (3.7) 



3.2. Quantum Liouville equation. We note that 

g a (t) := U(t, s)g a (s)U{s, t), t, s e R, (3.8) 

seen as a map from W 1,2 (M.) into W~ 1,2 (S.) is differentiable and solves the quan- 
tum Liouville equation (|1.25|) satisfying the initial condition g a (t)\ t = s = g a (s), 
provided g a (s) leaves W 1,2 (R) invariant. Indeed, using (|3.6p and (|3.7|) we find 

d 

— (g a (s)U(s,t)x,U(s,t)y) = 

i(U( S ,t)H a (t)x,g a ( S )U( S ,t)y)~i((g Q ( S )U( S ,t)x,U(s,t)H a (t)y) = 
i(H a (t)x, g a (t)y) - i(g a (t)x, H a (t)y), 

x,y e W /1 ' 2 (R), which yields 

d 

igl(Q a (t)x,y) = (g a (t)x,H a (t)y) - (H a (t)x, g a (t)y), 
x,y 6 W 1,2 (R), t.sel, 



5.5. Time dependent scattering. We set J7(t) := £/(£, 0), tsl and consider the 
wave operators 

:= s- lim U{t)*e- ltHo 

t — > — oo 

and 

Q+ := s- lim U(t)*e- ltH . 

t— * + OC 
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Proposition 31 Let Hp and H a (t), t £ R, a > 0, be given by (JT78]) - (jT7T6j) and 
(fl~2T|) - (fi~22|) . respectively. Then the wave operator I2_ and the limit 

R- := s- Km U{t)*{H D +T)- l U{t) (3.9) 

t — > — oo 

exist. Moreover, 

Ran ± (i7_) = Ker(i2_). (3.10) 



Proof. We start with (|3.9p . Let us introduce the time-dependent identification 
operator 

J D (t) -.^(H^ + Tr^HD + T)- 1 , teR. 

We have 

±U{t)*J D {t)e- UH °f= (3.11) 
iU(t)* {{H D + r)- 1 - [H a {t) + t)- 1 ) e~ itH °f + U {t)* j D {t)e~ itH f 
for / £ Sj where Jrj := 4rJo{t)- Hence we get 

U{t)J D {t)e~ UH °f - U( S yj D (s)e- lsH °f = (3.12) 
i j U{s)* {{H D + r)- 1 - (H a {r) + r)- 1 ) e~ ir H ° fdr + 

J S 

t 

irHr 



U(r)*J D (r)e- trH °fd S . 

Using (f3~4"| we find 

(H a {t) + t)- 1 = {H + t)- x I 2 Q b {H + r)-!/ 2 + (3.13) 
e at {H + T)~ l ' 2 Qi{e at + B^Q^H + t)- 1 ' 2 

where Qb is the orthogonal projection onto the subspace Ker(S). Note that Qg 
has rank 2. Taking into account (|1.23|) we get the representation 

(H D + t)- 1 = {H + t)- 1 ' 2 Q b {H + t)- 1 ' 2 . (3.14) 

By (l3~T3]l and ([3~T4]) we obtain 

{H D +r)- 1 -(H a (t)+T)- 1 = (3.15) 
-e at {H + T )-^ 2 Qi{e at + B)- l Q^{H + t)^ 2 . 

Since B is positive and invertible on Ker ± {B) we get the estimate 

\\{H d +t)- 1 -{H a {t)+T)- 1 \\ <e at \\QiB- l Qi\\, t £ R, a > 0. (3.16) 

Using again (|3.13[) we have 

j D {t) = ae at (H + T)- x l 2 Q^{e at + B)~ 2 BQb(H + r)-^ 2 {H D + t)" 1 . (3.17) 
This gives the estimate 

<ae at || Q^Qi || . (3.18) 
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Using (|3.12p . (|3.16|) and (|3.18p we prove the existence of the limit 



-itHr 



f2- := s- lim U(t)*J D (t)e 

t — > — oo 

In fact, the convergence is in operator norm: 

lim \\fi_ - U(t)*J D (t)e- ltHD \\ = 0. (3.19) 

t — > — oo 

Using the identity 

U(t)*J D (t)e- UHD - U(t)*e- ltH °(H D +T)- 2 
= U(t)*((H a (t) + r)- 1 - (H D + T)- x )e- itHD {H D + r)" 1 

and (|3.16|) we get the estimate 

\\U(tyj D (t)e- UH ° -U(t)*e- UH °{H D + T )- 2 \\ < e^WQ^B^Qil (3.20) 

which yields 

lim || /L - U(t)*e- ltHD (H D +t)- 2 \\ = 0, (3.21) 

t — > — oo 

for all t S K, a > 0. Since the wave operator fl_ exists we get the existence of 
s- lim U(t)*e- ltHD (H D +T)- 2 . 

t — > — oo 

Using that Ran((ffc> + r) -2 ) is dense in $3, we prove the existence of In 
particular, this proves that Q_ is isometric, i.e Q*_(2_ = I. 

Now let us prove that the operator in (|3.9p exists. Note that the norm con- 
vergence in p,19p yields the same property for adjoints: 

lim ||fr - e ltHo J D (t)*U(t)\\ = 0. 

t— oo 

In particular 

Q*_ = s- lim e ltH ° J D (t)*U(t). 

t — ► — OO 

In the quadratic form sense we get 

jU(t)*(H a (t) + r)- l U{t)f = u(ty lj t (H a (t) + r)- 1 j U(t)f 

f € S3, t e E, a > 0. Hence 

U(ty(H a (t) + r)- l U(t)f - U(s)*(H a (t) + T)~ 1 U(s)f 

drU{r)* i-^(H a (r) + r)- 1 Jf/(r)/ 
/ e S3, t, s e R, a > 0. By (33) we get 

4 (#*(*) + = ae- at (H + t)- x / 2 {I + e - at B)- 2 B(H + r)' 1 / 2 
at 
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which gives the estimate 

d_ 

dt 



{H a {t)+r)- 1 



<ae at \\Q^B-'Q^ 



(3.22) 



Hence R- exists, and we even have convergence in operator norm: 
lim - U{t)*(H a (t) +T)- 1 U{t)*\\ = 0. 

t — > — oo 

Taking into account the estimate (I3.16| we find 



lim \\R_ - U(t)*(H D + r)- 1 U(t)*\\ = 0. 

t—* — oc 



(3.23) 



In particular we have 



lim \\R 2 _ - U(tY(H a (t)+T)- 2 U(t)*\\ 

= lim \\R 2 _ - U(t)*(H D +T)- 2 U(t)*\\ = 0. 

t — > — OO 

Using the identity 

J D (t)* = ((H D + t)- 1 - (H a (t) + t)- 1 ) (H a (t) + t)- 1 + (H a (t) + T )- 
and taking into account the estimate (|3.16p we obtain 



lim \\R±-U(t)*J D (t)*U(t)\\=a. 

t — > — oo 



Hence we find 



f2*_ = s- lim e ttHD J D (t)*U(t) 

t — * — OO 

= s- lim e ltHD U(t)U(t)*J D (t)*U(t)= s- lim e ltHD U(t)R 2 _ 

t — > — oo t — > — OO 

Mr 



(3.24) 



which shows in particular that the limit lim t _ > _ 00 e D U(t)f exist for elements 
/ £ Ran(i?_). More precisely: 



lim e itHD U(t)R_f = S1*_R_f 

t— > — OQ 



(3.25) 



for all /. 



We are now ready to prove (|3.10p . Assume that / _L Ran (.!?_). Then using 
the definitions, the unitarity of U(t)*, and (|3.23[) we obtain: 

0= (/, fJ_(H D + T)- 1 g) 
= lim (U{t)*{H D + r^U^f, U{t)*e~ ltH °g) = (J2_/, Q.g) 

t — > — OO 

for g e f). Hence / 1 R an(f2_) implies R-f _L Ran(/2_) = Ker(J2*). Thus 
Q*_R-f = 0. Using (I3T251) : 

0= lim \\e itH °U(t)R.f\\ = \\R_f\\, 

t — ► — OO 
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thus / € Ker(i?_ ) . We have thus shown that Ran ) C Kcr(i?_ ) . Conversely, 
choose / £ Ker(i?_). We have (use lpT23]l ): 

(/, f2-(H D + r)-^) = lim (/, U(tye-* tH ° (H D + r^g) 

I— ' — CO 

= lim (U(ty(H D +T)- 1 U(t)f,U(tye- itH °g) = (R-f,f2_g) = 0, 

t — >— oo 

for all g. Thus i?* / is orthogonal on a dense set (domain of Hr>), thus equals 

zero. Therefore Ker(i?_) C Ran ± (J7_) and ()3.10|) is proved. 

Remark. Note that J7_ would be unitary if one could prove that Ker(i?_) = 0. 

Proposition 32 Let H and H a (t), t e M., a > 0, be given by (|l.l|) -(|1.4 |l and 
(jl.21[) - (|1.22[) . respectively. Then the wave operator f2 + exists and is unitary. 

Proof. We introduce the identification operator 

J(t) (^(tJ + r^^ + Tr 1 . tern.. 

In the quadratic form sense we get that 

±U(tyj{t)e- aH f = (3.26) 
iU(ty((H + t)- 1 - (H a (t) + T )- x )e* H f + U(tyj(t)e- itH f, 
iel, where 

j(t) := |j( t ). 
Taking into account (|3.4p we find 

(F+r)- 1 -( J ff Q (t)+r)- 1 =e- at (H+T)- 1/2 (I+e- at B)- 1 B(H+T)- 1/2 , (3.27) 
i £ R. Hence we have the estimate 

IK// + T-)- 1 - <e- Qt ||£||, t€R. (3.28) 
Moreover, we get 

j(«) = |-(ff Q (t)+T)- 1 (F + T)- 1 (3.29) 
ar 

= ae- at (H + t)- 1 ' 2 ^ + e- at B)- 2 B(H + rT 3 / 2 
which yields the estimate 

m 



<ae~ Qt ||£||, te 



Hence the strong limit 



Q+ : = s- lim U(tyj(t)e 



-itH 



exists. Moreover, the convergence is also true in operator norm: 

lim 0+ - U(t)*J(t)e- itH \\ = 0. (3.30) 

t — » + oo 
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Using the identity 

U{t)*J{t)e- ltH - U{t)*e- UH {H + t)~ 2 = (3.31) 
U(t)* {(H a (t) + r)- 1 -(H + r)- 1 ) e~ ltH (H + r)" 1 

and taking into account the estimate ()3.28j) we obtain 

lim \\f2 + -U(t)*e- UH (H + T)- 2 \\=0. (3.32) 

t — >+oc 

Hence /2+ exists on a dense domain and is isometric, i.e. Q+fi + = I. 

Let us now prove that .!?+ is unitary. Since (|3.30p holds also true for adjoints 
we find 

lim \\ni-e itH J(t)U{t)\\=0. (3.33) 
Hence, we have the representation 

H* + = s- lim e ltH J(t)U(t). (3.34) 

Furthermore, in the quadratic form sense we have 

j t U{t)*{H a {t) + r)- l U{t) = U(ty [j t {H a {t) + r)- l] jU(t) = 
ae- at U(t)*(H + t)- 1 / 2 ^ + e - at B)- 2 B(H + T)- 1 ' 2 U{t). 
Hence we get 

U(t)*(H a (t) + T)- x U{t) = (H a (0) + t)- 1 (3.35) 

+a f e- at U(t)*(H + t)- x ' 2 {I + e - at B)- 2 B(H + T)-^ 2 U{t)dt. 
Jo 

Using the estimate 

U(t)*(H + t)- 1 / 2 (I + e- at B)- 2 B(H + T)-^ 2 U{t) < \\B\\, t e R, 
we find that the following weak integral exists and defines a bounded operator: 

/>oo 

a dte- at U{t)*{H + t)- 1 / 2 (I + e- at B)- 2 B{H + T )- 1 ' 2 U{t). (3.36) 



Moreover, by the Cook argument it also implies the existence of the limit 
R + := s- lim U{t)*(H a {t) + r) -1 ?/(*). 

t — >+oo 

In fact, the convergence takes place in operator norm: 

lim \\R+ - U(t)*(H a {t) +TT 1 [/(i)|| =0. 

t — >+oo 

Taking into account the estimate (|3.28p we obtain 

lim \\R+-U{t)*{H + T)- 1 U{t)\\=0. 

' + 00 
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which yields 

lim \\Rl-U(t)*J(t)U(t)\\ = 0. 

By (pT33"l) we get 

^Hm^H^-e^C/^ii^ll =0 

which shows the existence of fi+f = s-lim t ^ +00 e ltH U(t)f for / s Ran(i? + ). 
Now in order to prove the unitarity of J? + it is enough to show that Ran(i? + ) 
is dense in $). Let us do that. 
From (|3.35j) we obtain 

R+ = {H a (0)+ T )- 1 + 

a dt e- at U{t)*{H + t)- 1 ' 2 ^ + e - at B)- 2 B(H + r)-^ 2 U(t), 
Jo 

which by the positivity of the integral it gives < (H a (0) + r) _1 < R + . Hence 
Ker^ 2 ) C Kcr((ff Q (0) + r)" 1 / 2 ), 

thus 

Ran^^y 2 ) C Ran ± ((if Q (0)+r)- 1 / 2 ) =0. 

Thus we get that Ran(i?V 2 ) is dense in fj which yields that Ran(i?^) is dense 
in f). Therefore we have the representation = s-lim t ^ +00 e UH U(t) which 
proves that Q + is unitary. 

3.4- Time- dependent density operator. Now we are ready to write down a solu- 
tion to our Liouville equation (11.25|) which also obeys the initial condition at 
t = — oo. Let us introduce the notation: 

g a {0) := fl-Q D Q*_ 

which defines a non- negative self-adjoint operator. Here go is given by (|1.17[) - 
(|1.19[) . In accordance with (|3.8p . the time evolution of g a (0) is given by 

g a {t) = U{t)g a {Q)U{t)* = U{t)Q^Q D Q*_U{t)* 1 iel, (3.37) 

where we have used the notation U(t) := U(t,0) and the relation U(0, t) = U(t)*, 
tgl. We now show that the initial condition is fulfilled. 

Proposition 33 Let Ho and H a (t), t G R, a > 0, be given by (|1.8p - (|1.16[) and 
(|1.21j) - ([1.22|) . respectively. If go is a steady state for the system {Sj,Ho} such 
that the operator go '■= {Ho + t) Qd is bounded, then 

lim \\g D -g a (t)\\ =0. (3.38) 
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Proof. We write the identity: 

U{t)Q- QD n*_U{t)* = U{t)Q^{H D + T)- 2 g D {H D +T)- 2 Q*_U{t)\ (3.39) 
t e M. Taking into account (|3~2T1) we find 

u(t){2-g D ntu(t)* = u{t)f2-g D s}*_u(ty. (3.40) 

From (|3.19[) we get 

lim \\U(t)h-Q D h*_U{t)* - J D (t)e- UHD Q D e ltHD J D {ty\\ (3.41) 

t— * — OO 



lim ||J7(i)/2_gD/2* J7(i)* - Jr. •/£>(*)! = 0. 

£ — > — OG 



Using (|3.16|) we get 



lim ||JD(t)?DJD(*)*-(#Z5+Tr 2 ?D(#D+T)- 2 || (3.42) 

£— *■— 00 

= lim || J D (t)gD J D (t)* - toll =0 

t—> — OO 



Taking into account (|3739|) - (j3T42|) we prove (|3T38|) . 



5.5. Large time behavior on the space of absolute continuity. We now are ready 
to prove the result announced in (|1.27p . 

Proposition 34 Let H and H a (t), t 6 R, a > 0, fee giwen 6y (|1.1[) - (|1.4[) and 

(fL2T1) - (fl~22|) . respectively. Let W- be the incoming wave operator as defined in 
(|1.26p . If qd is a steady state for the system {ft^Hf)} such that the operator 
qd '■= {Hd +t) 4 (?Z5 is bounded, then 

s- lim g a {t)P ac (H) = W-Q D W\ (3.43) 

t — >+oo 

Proof. Let us assume that the following three technical results hold true: 

s- lim (U(tY - e ttH )P ac (H) = 0, (3.44) 

(H d +t)- 2 {[}*_- 1) is compact, (3.45) 

and 

s- lim (H D + t)- 2 {Q*_ - I)e ltH P ac (H) = 0. (3.46) 

t — >+oc 

We will first use these estimates in order to prove the proposition, and then 
we will give their own proof. 

We write the identity: 

U(t)g a (t)U(tyP ac (H) 
= U{t)f2_{H D + T)- 2 g D {H D + T )- 2 f2*_U{tyP ac (H) 
= U(t)n-(H D +t)- 2 q d (H d +t)- 2 [2*_ (U(t)* - e itH )P ac (H) 
+U(t)Q-(H D + t)- 2 q d (H d + r)- 2 (rr - I)e ltH P ac (H) 
+U(t)n_(H D + r)- 2 e ltHD Q D (H D + T)- 2 e~' ltHD e ltH P ac (H). 
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Taking into account ()3.44() - (|3.46|) . and using the completeness of W~ which yields 
Wl = s-lim^-oo e ttHD e- ltH P ac (H), we get: 

s- lim U(t)g a (t)U(t)*P ac (H) 

t — > + oo 

= s- lim U(t)Q_(H D + T)- 2 e ltHD g D (H D + t)~ 2 W1. 

t — >+oo 

Since — I) (Hp + t)~ 2 is also compact (its adjoint is compact, see (|3.45p ). 
we have: 

1-2 AtHr 



s- lim (!?_ - I)(H D + T)- z e ltMD P ac {H D ) = 0. 



Thus: 



s- lim U(t)g a (t)U(t)*P ac (H) = 

t—>-\-oo 

s- lim U{t)e ltHo (H D + t)- 2 q d {H d + t)- 2 W1 

= s- lim U(t)e ltHD g D Wl = s- lim U {t)e ltH P ac {H)W„ Q D W*_. 

t — >+oo t — > + oo 

Finally, we apply (|3.44|) once again, and (I3.43|) is proved. 

Now let us prove the three technical results announced in (|3.44p - (|3.46p . We 
start with (j3~44) . 

We have the identity: 

(U{t)* - e UH ){H + t)- 2 = (U(t)*(H + T)- 2 e- ltH - (H + r)- 2 )e ltH . (3.47) 

Then by adding and subtracting several terms we can write another identity: 

(U(t)* -e UH )(H + t)- 2 
= [u(t)*(H + r)- 2 e- ltH - e ltH (3.48) 

+ J(0)} e ltH g + { J(0) -(H + r)~ 2 } e ltH . (3.49) 

By 4222) we get 

lim \\U(t)*(H + r)- 2 e- ltH - f2+\\ = 0. (3.50) 

which shows that (|3.48[) tends to zero as t — » +oo. Next, from (13. 4| . (13.26|) and 
dHUD we get 

U{t)*J(t)e~ ltH - J(0) 

= t / ds e- as U(s)*(H + T y 1/2 (I + e- at B)- 1 B(H + r)- 1/2 e- lsH 
Jo 

+a [ ds e- as U(s)*(H + t)-^ 2 (I + e- as B)- 2 B(H + T)- 3 / 2 e- lsH , 
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which yields 
i?+ - J(0) 

/"OC 

= i / ds e- as U(s)*(H + t)- 1 ' 2 ^ + e~ as B)~ 1 B(H + T)- 1 l 2 e~ lsH 
Jo 

+a f ds e- as U(s)*{H + t)- 1/2 (I + er as B)- 2 B(H + r)^^ e~ lsH . 
Jo 

Since B is a compact (rank 2) operator, we get that i? + — J(0) is a compact 
operator. This fact immediately implies (via the RAGE theorem): 

s- lim (J7+ - J(0))e itH P ac (H) = 0. (3.51) 

t^>-\-oo 

Furthermore, we have the identity: 
J(0)-(H + r)- 2 
= ((H a (0) + T)- 1 -(H + r)- 1 )(H + r)- 1 

= -(H + t)- 1 / 2 ^ + B)~ 1 B(H + r)~ 3 / 2 , 
which gives that J(0) — (H + t)~ 2 is compact. Thus: 

s- lim (J(0) - [H + T)- 2 )e aH P ac {H) = 0. (3.52) 

Taking into account (|3"30j) , (|33T|) and ([3321) we find 

s- lim (U(t)* - e ltH )(H + r)- 2 P ac (H) = 

t — > + oo 

which proves (|3.44|1 . 

Next we prove (|3.45[) (the estimate p.46[) is just an easy consequence of (|3.45[) 
via th e RA GE th eorem). Using (pOT]) we have f2_ {H D +I)- 2 = J?_ . From ((21121), 
([345]) and ((347)) we obtain: 

J D (0)-U(t)*J D (t)e- itHD 

r° 

= -i dse as U{s)*(H + T )- 1/2 Q^{e as + B^Q^H + r)- 1/2 e- lsHD 

r° 

+a J ds e as U(s)*(H + r)- 1/2 (e as + B)- 2 B(H + t)- 1 I 2 (H d + r^e'" 110 

and together with the fact that Q% is a rank 2 operator we find that i?_ — Jd (0) 
is compact. Hence f2*_ — Jd{0)* is compact, too. 
Moreover, using (|3.15|) we get 

Jd(0)-(H d +t)- 2 = 

((H a (0) + r)- 1 - (H D + t)- x ){H d + r)- 1 = 

-e at (H + T)- l ' 2 Q^(e at + B^Q^H + r)-^ 2 (H D + r)" 1 

which shows that Jd(0) ~ {Hd + t)~ 2 is compact. Hence Jd(0)* — {Hd + t)~ 2 
is compact. 

Now use the identity: 

(H D + t)- 2 {[2*_ -I) = {Q*_ - J D (0)) + (J D (0) - (H D + r)" 2 ), 

which proves (|3.45[) . Finally, (|3.46p follows from (|3.45p and the RAGE theorem. 
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3.6. The main result. We are now ready to rigorously formulate and prove our 
main result, announced in the introduction: 

Theorem 35 Let H andH a (i), t e R, a > 0, be given by ([L~T])-([0]) and (fl~2"Tjl - 
(|1.22|) , respectively. Let W~ be the incoming wave operator from (|1.26p . Further, 
let Eh(-) and {Xj}^^ be the spectral measure and the eigenvalues of H . If qd is 
a steady state for the system {Sj, Hp} such that the operator qd := (Hp +r) 4 qd 
is bounded, then the limit 

1 f T 

g a := s- hm — / dtg a (t) = (3.53) 

N 

^EaUXj^ScQDS^EH^Xj}) + W-q d W1 

3 = 1 

exists and defines a steady state for the system {.f), H} where S a := 

Remark 36 We stress once again that only the part corresponding to the pure 
point spectrum g^, := Ylj=i ^H({Aj})S a QDS^EH({\j}) of the steady state g a 
depends on the parameter a > while the absolutely continuous part g^ c ■= 
W-goWt does not. Note that with respect to the decomposition Sj = F) P (H) © 
Sj ac (H), one has Q a = Q p a ® g a a c . 

Proof. By Proposition [M] we have 

1 f T 

s- lim - / dt g a {t)P ac {H) =W-g D W*_. (3.54) 



1 ' T 



In particular, this yields: 



s- lim - / dt P s {H)g a {t)P ac {H) = 0, 

where P S (H) is the projection onto the singular subspace of H. Now we are 
going to prove 

s- lim i / P ac {H)g a (t)P s {H)dt = 0. (3.55) 



T 



By (j3~32|) we find 



lim \\U(t)f2_g D n*_U(t)*(I + H)- 2 e- ltH - U(t)n_g D f2*_n+\\ =0, 

t — >oo 

which yields 

lim \\U(t)Q-Q D n*_ U{t)*{I + H)- 2 - U{t)f!-g D n*_n + e ltH \\ = 0. 

t — >oo 

Let Xj be an eigenvalue of H with corresponding to an eigenfunction <j)j. Then 
lim \\U{t)Q_g D Q*_ U{t)*{I + H)- 2 4>,j - e ltx >U{t)Q_g D [2*_ fl+faW =0. 

t— *oo 
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This and the unitarity of f2\ give: 

\-2 



lim \\g a (t)(I + H)- 2 <t>j 



e it\ ie -itH 



= o 



or 



lim \\g a (t)^j 



n* + n_Q D n*_n + (j) :j \\ = o. 



(3.56) 



Hence we have 



dt P ac (H)g a (t)cj)j = 



(i 



T 



dt 



P ac (H)f2lf2_g D f2*_f2+(j) J . 



We use the decomposition 



I J P ac {H)g a {t)4>, 

i r T 

- dt e ttx i e - ltH E H (\\ -Aj\< e)P ac (H)nin_g D f]*_a 
T Jo 



T 
1 
T 



dt 



itXj -itH 



E H (\X - Xj\ > e)P ac {H)Q* + Q^g D Q*_Q+(, 



If e is small enough, then Eh(\X — Xj\ < e)P ac (H) = 0. This yields the estimate: 



P ac (H)g a {t)<t> 3 



1 




T 


L 




2 


< 






T 



which immediately shows that 



s- lim i / dt P ac {H)g a {t)ct) 3 ; = 0, 

T->oo i J Q 



(3.57) 



and (|3.55p is proved. Next, from (|3.56p we easily obtain: 



s- lim — 



dt P s {H)g a {t)cj) = E H ({Xj}) S a g D S*E H ({Xj}). (3.58) 



Now put together (|3"3I)) . (|3"35j) . (j3~57) and (|3351) . and the proof of (|3~53l> is 
over. Now the operator g a is non- negative, bounded, and commutes with H. 
Hence g a is a steady state for {Sj, H}. 

Corollary 37 Let H and H a (t), t € K, a > 0, be given by (pTT]) - (pT4 jl 

and (|1.21[) - (|1.22[) . respectively. Then with respect to the spectral representation 
{L 2 (R, f)(A), ^), M} of H the distribution function {/5 q (A)}a£R of the steady state 
g a is given by 



Pa(X) 



f&(A - /i fc ), 
ffbiX-fib) 

I fa(A-AO 



A eM\ff(ff) 
A = Aj, j = l, 
A 6 K,w a ) 

A e [u a , oo) 



(3.59) 



w/iere p aj - := (S a <f>j,^j), j = 1, 2, . . . , JV. 
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Proof. Using the generalized Fourier transform (|2.2|) - (|2.4p one has to consider the 
operator <P~ 1 g a <P : L 2 (R, f)(A), v) — > L 2 (R, f)(A), v). Using the representations 

'pQa'p- 1 = 'Pg 1 ^ 1 + <Pg a a c <P' 1 

and 

'Pg^'P^ 1 = 'PW-g'gWl'p- 1 = $W-$ r - 1 &Q'%&- 1 &Wl$- 1 , 
Qd ■= Qa © Qb, we get 

M(p° c ) = $W-&- 1 M(p'£)iI'Wl<p- 1 . 

where 



0, XeR\<r ac (H) 

A e [vb,v a ) 

A e [v a , oo). 



. fb(A — Ub) A G \vh, v a ) 

PdW--= / \ n \ n \ (3.60) 
ft (A - /i 6 ) > 



fo(A-// a 
Taking into account (|2.5p we prove (|3.59p . 



3. 7. The case of sudden coupling. Let us compare our results with following 
model Assume that our system is not coupled for t < and suddenly at 
t = the system becomes fully coupled. In a more mathematical manner this 
can be modeled by the following family of self-adjoint operators: 



H(t) := 

To the evolution equation 



H D t<0 
H f > 0. 



d 

i— u (t) = H(t)u(t), t e M, 

it corresponds a unique unitary solution operator or propagator 
{U(t, s)}( M)eR><R given by 

U(t,s) := U(t)U{s)-\ t,seR, 

where 

U(t) 

The time evolution of the density operator is given by 
Qoo (t) := U(t) QD U(t)*, t > 0. 
Clearly, lim^-oo ||^oo(t) — toll = 0. Then using the identity: 
g(t) = e- ltH e* tHD g D e- ltHD e ltH , t > 0, 



e- lfH , t > 0. 
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we immediately get that 



Hence we find 



s- 

T 



s- lim goo {t)P ac (H) = W- BD W*_. 

t— »+oo 



lim 1 / dt Qoo (t)P ac (H) = W.QdWI 
^+00 1 J 



As above, we can show that 

rT 



s- lim 1/ dtP ac (H) Qoo (t)P s (H) = 



and 

N 



S- 

T 



lim 1/ ^P s ( J ff)g 00 (t)P s ( J ff) = ^£; ff ({A J }) to ^({A,}). 
-+°° 1 Jo ~i 



Hence we find 

N 



lim ^ f dt Qoo {t) = Y / E H {{\ j }) QD E H ({\ j }) + W-QdWI 



4. The stationary current, the Landau-Lifschitz and the 
Landauer-Biittiker formula 

There are by now several proofs of the Landauer-Biittiker formula in the NESS 
approach (see [4"I|27|). and in the finite volume regularization approach (see [T51 
H4"1H"2"] ). Here we give yet another proof in the NESS approach. In fact, we will 
only justify the so-called Landau-Lifschitz current density formula (see (|4. 14[) in 
what follows), which was the starting point in [6 for the proof of the Landauer- 
Biittiker formula (see Example 5.11 in that paper). 

Let us start by defining the stationary current, in the manner introduced in 
[4]. Let -q > 0, and choose an integer N > 2. Denote by \b the characteristic 
function of the interval (b, 00) (the right reservoir). Without loss of generality, 
let us assume that H > 0. 

Definition 41 The trace class operator 

3 {f 1 ):=i[H{l + i 1 H)- N :X b] (4.1) 

is called the regularized current operator. The stationary current coming out of 
the right reservoir is defined to be 



3 a := lim Tr (g a j(rj)). 



(4.2) 
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Now a few comments. The current operator is trace class because we can write 
it as 

i[H{l + r)H)~ N - H D (1 + r,H D )- N , X b] 

which clearly is trace class. Then since g v a does not contribute to the trace in the 
definition of the current, we will focus on the clearly a independent quantity: 

3 = lim Tr {W- goWlP^^j^)} . (4.3) 

We start with a technical result: 

Lemma 42 Let x a bounded, compactly supported function. Then the operator 
x(l + H)~ 2 is trace class. 

Proof. Choose a smooth and compactly supported function x such that XX — X- 
Write 

x(i+h)- 2 = xx(i+h)- 2 = xii+Hy'xii+Hy'+xii+Hy^H^Ki+H)- 1 . 

Since x is smooth an compactly supported, the operator xli^+H)^ 1 
is Hilbert-Schmidt. The operators x(l + H) -1 and x(l + are also Hilbert- 

Schmidt, thus x(l+-ff) _2 can be written as a sum of products of Hilbert-Schmidt 
operators, therefore it is trace class. 

Next, let us now prove that we can replace the sharp characteristic function 
in the definition of the current with a smooth one. Let c > 6+1. Choose any 
function <j) c G C°° (E) such that 

O<0 C <1, <t>c{x) = 1 if x > c+1, suppO c ) c (c- l,oo). (4.4) 
Then let us prove the following identity: 
Lemma 43 

Tr {W-QdW! P ac (H)j(r])} = (4.5) 
zTr {W-g D WlP ac (H)[H(l + V H)- N , C ]} . 

Proof. First, the commutator [H(l + rjH)~ N , <j> c ] defines a trace class operator; 
that is because now [H,cf> c ] = + ^|). and i 1 + ^) _1 [^0c](l + 

rjH)^ 1 is a trace class operator (we can write it as a sum of products of two 
Hilbert-Schmidt operators). We also use the identity 

W-g D WlP ac {H) = W-g D (l + H D )W*_P ac (H){l + H)- 1 (4.6) 

which is an easy consequence of the intertwining property of W— ■ 
Second, (|4.5|) would be implied by: 

Tr{W-Q D WlP ac (H)[H(l + r}H)- N ,(l> c -Xb}} = 0- (4.7) 

We see that <p c — Xb has compact support. If we write the commutator as the 
difference of two terms, both of them will be trace class. The first one is 

W-g D WlP ac (H)H(l + V H)- N (cj )c - X b) 

= W-g D (l + H D ) 2 WlP ac {H)H{\ + vH)- N {(l + H)- 2 (cf> c - X b)} 
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and the second one is 

W-Q D WlP ac {H){cjy c - Xb )H{\ + V H)- N 

= W-g D (l + H D ) 2 W*_P ac {H){{\ + H)- 2 (^ c - X b)}H(l + V H)- N . 

Now according to Lemma |4"2I (1 + H)~ 2 {(f) c — Xb) is a trace class operator. Thus 
the two traces will be equal due to the cyclicity property and the fact that H 
commutes with the steady state. 

We can now take the limit rj \ 0: 
Lemma 44 The operator (1 + H)~ 2 [H, <f> c ] is trace class, and 

3 = iTr{W-Q D WlP ac (H)[H,0 c }} , (4.8) 

independent of </> c . 

Proof. Note that [H,<t> c ] = -^(Z-^ftc ~ 4>"), where both <j>' c and <#.' are 
compactly supported. Using the method of Lemma 22] one can prove that 
(1 + H)^ 2 -^(/)' c is trace class, hence (1 + H)~ 2 [H, cj) c ] is trace class. Thus 
W-QDWlP ac (H)[H, cf) c ] is trace class since we can write 

W- eD WlP ac {H)[H^ c ] = W-g D (l + H D ) 2 WlP ac (H)(l + H)- 2 [H,4> C }. 

In fact, using trace cyclicity one can prove that 

Tv{W^g D W:P ac (H)[H,cb c ]} = (4.9) 
Tr {W-g D (l + H D ) 2 W*_P ac {H)(l + H)- X [H, cf> c ](l + H)- 1 } 
= -TY{W-Q D {\ + H D ) 2 WlP ac {H)[(\ + H)-\cp c ]} . 

This last identity indicates the strategy of the proof. Write: 

Tr {W- QD W*_P ac {H)[H(l + nH)- N ,cb c }} (4.10) 
= Tr {W-g D (l + H D fWlP ac {H){l + H)- 2 [H(l + T]H)- N , c ](1 + H)- 1 } . 

Now it is not so complicated to prove that (1 + H)~ 2 [H(1 + t]H)~ n , </> c ](l + -ff) _1 
converges in the trace norm to (1 + H)~ 2 [H, </> c ](l + H) -1 when rj \ 0; we do 
not give details. Now use (|4.5p and take the limit; we obtain: 

3 = iTr {W^g D (l + H D fW*_P ac {H)(l + H)- 2 [H, C ](1 + H)- 1 } 
= iTr{W^Q D WlP ac (H)[H, ( j) c ]}, (4.11) 

where in the last line we used trace cyclicity. 

The remaining thing is to compute the trace in (|4.8p using the spec- 
tral representation of H. Let us compute the integral kernel of A := 
iW.g D WlP ac {H)^ (-£<t>' c - (f>' c £) in this representation. We use ([3^0]) . 
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where we denote the diagonal elements of £>£?(A) by g D c (X) pp (the other entries 
are zero). We obtain: 

A(X,p;X',p') = (4.12) 

-^S c (A) PP ^0 P (x,A) f^tf>' e {x) + tf>' c (x)-^j $ p ,(x,X')dx 

= -7^9dWpp I 4>' c {x){Mx,X)^ p ,{x,\') - 4>' p (x, \)<t> p r(x,A')}dx, 

where in the second line we integrated by parts (remember that <fi' c is compactly 
supported). 

In order to compute the trace, we put A = X',p = p', and integrate/sum over 
the variables. We obtain: 

3= / K(x)j(x)dx, (4.13) 

JR 

where 

1 poo 

j(x):=— ^a c (A) pp Im{^(x,A)0 p (z,A)}dA (4.14) 

mb Jvt p 

is the current density, which can be shown to be independent of x (the above 
imaginary part is a Wronskian of two solutions of a Schrodinger equation, see 
[6] for details). But J R 4>' c (x)dx = 1 for our class of cut-off functions, therefore 
the stationary current equals the (constant) value of its density. The Landauer- 
Biittiker formula follows from the Landau-Lifschitz formula (|4.14[) as proved in 

®- 
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